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Abstract 

We consider a class of differential equations, x + ~/x + g(x) = /(w£), with u> £ WL d , describing 
one-dimensional dissipative systems subject to a periodic or quasi-periodic (Diophantine) forcing. We 
study existence and properties of the limit cycle described by the trajectory with the same quasi- 
periodicity as the forcing. For g(x) — x 2p+1 , p £ IN, we show that, when the dissipation coefficient 
is large enough, there is only one limit cycle and that it is a global attractor. In the case of other 
forces, including g(x) = x 2p (with p = 1 describing the varactor equation), we find estimates for the 
basin of attraction of the limit cycle. 

1 Introduction 

Consider the ordinary differential equation 

x + j± + x 2p+1 = f(ut), (1.1) 
where p £ M, u> £ M, d is the frequency vector, /(i/>) is an analytic quasi-periodic function, 

with average (/) = fo ^ 0, and 7 > is a real parameter (dissipation coefficient). Here and henceforth 
we denote with • the scalar product in TH d . By the analyticity assumption on / there are two strictly 
positive constants F and £ such that one has \f v \ < Fe - ^'"' for all v £ r L d . 

If d > 1 we shall assume a Diophantine condition on the frequency vector u>, that is 

|w -v\>C Q \v\- T Vi/eZ d \{0}, (1.3) 

where \v\ = \v\\ = + . . . + \vd\, and Co and r are positive constants, with r > d — 1 and Co small 
enough. Note that for d = 1 the condition (|1.3|) is automatically satisfied for all us. 

In this paper we want to show that for 7 large enough the system (|l.lf) admits a global attractor 
which is a quasi-periodic solution with the same frequency vector u> as the forcing /. This will be done 
in two steps: first we prove that for 7 large enough there is a quasi-periodic solution Xo(t) with frequency 
vector u) (cf. Theorem^! 11 Section^, second we prove that, again for 7 large enough, any trajectory is 
attracted by Xo(t) (cf. Theorem |2 in Section |3J|. 

In particular, this solves for the system i|l.l|) a problem left as open in [7|. Indeed in [7j we considered 
a class of ordinary differential equations, including (|1.1|) . and proved existence of a quasi-periodic solution 
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with the same quasi-periodicity as the forcing, but we couldn't conclude that this was the only solution 
with such a property. The result stated above gives an affirmative answer to this problem for the 
system (f 1 . 1 1) , by showing that the quasi-periodic solution xo(t) is unique; cf. Theorem in Section 0] 

This uniqueness result holds for the more general systems studied in \7\ , including the resistor- inductor- 
varactor circuit, or simply varactor equation, described in In that case the solution Xo(t) is not a 

global attractor any more, but it turns out to be the only attractor in a neighbourhood of the solution 
itself. 

More precisely the situation is as follows. We can consider systems described by 

x + jx + g{x) = f(ut), (1.4) 

where / is given by 1|1.2|) and g is an analytic function. Studying the behaviour of the system (11.411 
for 7 large enough suggests to introduce a new parameter e = I/7, in terms of which the differential 
equation (|1.4(l becomes 

ex + x + eg(x) — ef(ujt), (1.5) 

and study what happens for e small enough. 

If we assume that there exists Co £ E such that g(co) = fo and g'(co) := d x g(co) ^ 0, then the 
system (|1.5H admits a quasi-periodic solution xo(t), analytic in t, with the same frequency vector u> as 
the forcing /, and furthermore xo(t) = cq + 0{e). This was proved in 0, where the solution xo(t) was 
explicitly constructed through a suitable summation of the perturbation series 

00 

=^2e k x {k) (ujt), xWty) = e ""M fc \ (1-6) 

fc=0 I/£Z J 

for a function formally solving the equations of motion. 

As a drawback of the construction we were not able to prove any uniqueness result about xo(t). In fact, 
in principle, there could be other quasi-periodic solutions near xo(t), possibly with the same frequency 
vector u}. Neither we could exclude existence of other solutions reducing to Co as e — ► or even admitting 
the same formal expansion 1)1. 6f) in powers of e. In this papers we get rid of these possibilities, and we 
prove that there exists, in the plane (x,x), a neighborhood B of the point (co,0) where (xo(t), &o(t)) is 
the only stable solution of (|1.5I) . Moreover it turns out to be asymptotically stable, that is it attracts any 
trajectory starting in B. Therefore, this allows us to formulate a strengthened version of the theorem 
proved in [7j; cf. Theorem 0] in Section 0J 

In general the neighbourhood B can be very small. In specific cases one can look for improved 
estimates of B. In particular we can consider the case g{x) = x 2 , corresponding to the varactor equation 
studied in |3] . In that case we can give a good description of the basin of attraction of the quasi-periodic 
solution; cf. Section^ This provides a complementary result to the analysis performed in 0, where the 
set of initial data generating unbounded solutions was extensively studied. The analysis can be extended 
to any function x 2p , p£ I. We also construct a positively invariant set containing the attracting periodic 
orbit which allows us, together with the results of Section [3 to obtain a larger set estimating the basin 
of attraction (with area growing linearly in 7); cf. Theorem [3] in Section [H] 

Finally we show that in the case g{x) — x 2p , p £ M, there are unbounded solutions which blow up in 
finite time, and we discuss the implications for the varactor equation in |3J; cf. Theorem in Section 

More formal statements of the results will be formulated in the forthcoming sections. Some open 
problems will be discussed at the end. Here we confine ourselves to noting that, while in the case of 
periodic forcing standard techniques, like those based on Poincare sections could be applied, this 

is not the case for quasi-periodic forcing, where no Poincare maps can be introduced. 

The rest of the paper is organised as follows. Sections [3 and |31 are devoted to the global study of the 
system (|1.1|) for 7 large enough, whereas in Section 0] we draw the conclusions, and we pass to the study 
of the system (11.411 . In Section [3] we specify the construction envisaged for any g(x) in (|1.4|) to the case 
g{x) — x 2p , with particular care for p — 1, and we look for an estimate of the basin of attraction of the 
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attracting periodic solution. Then, in Section H3 we use the techniques of |3] to improve the estimate of 
the basin of attraction. In Section we show that, again in the case g(x) = x 2p , there are unbounded 
solutions which blow up in finite time. Finally in SectionElwe mention some open problems, and possible 
future directions of research. 



2 Existence of the quasi-periodic solution 

First we show that for 7 large enough there exists a quasi-periodic solution xo(t). 

Theorem 1 Consider the equation with f a non-zero average quasi-periodic function analytic in 

its argument and with ui satisfying the Diophantine condition fi.ffi . There exists 70 > such that for 
all 7 > 70 there is a quasi-periodic solution xo(t) with the same frequency vector as the forcing term. 
Such a solution extends to a function analytic in I/7 in a disc T> of the complex plane tangent to the 
imaginary axis at the origin and centered on the real axis. Furthermore, Xo(t) = a + 0(1/7), where 

a = / 1/(2p+1) /o. 

Proof. We can apply the results of Section 7 in [7|. If we set g(x) — x 2p+1 , then <?(co) = fo yields 
Co = /o P j so that g'(co) ^ as by assumption one has fo 7^ 0. Both the existence of the analyticity 
domain T> and the form of the solution itself follow from the analysis in 0. ■ 

If 7 is large enough, say 7 > 7 > 70, then the solution xo(t) is of definite sign. In the following we 
shall assume that this is the case: hence xo(t) ^ for all t g R. 



3 Convergence to the quasi-periodic solution 

Given the quasi-periodic solution xo(t) one can write x(t) = xo(t)+^(t), with satisfying the differential 
equation 

t + rt + ZF(Z,x o (t)) = 0, (3.1) 



where we have defined 



2p / \ 

F& x) := i ((x + 2p+1 z 2p+1 ) = E r\ + 1 ) e p+1 -^ J . (3.2) 



We can write (|3.1|) as 



^ y ' (3 3) 

V = -TV-ZF(t,xo{t)), 



that is z = ${z), if we define z = (£,y) and $(2) = (y, —jy — £_F(^,x (t)). We denote by tp(t, z ) the 
solution of l|3.3l) with initial datum z . Define also P(£,t) := F(£,Xo(t)) and 0(f) := F(l;,a) and set 
R(^t):=P(U)/Q(0- 

Here we prove the following result. 

Theorem 2 Consider the equation with xa(t) the quasi-periodic solution of given in Theorem 

1. There exists 71 > such that for all 7 > 71 all trajectories in phase space converge toward the origin 
as time goes to infinity. 

The proof will pass through several lemmata. 
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Lemma 1 Assume 7 > 7 so that xn(t) exists and xo(t) 7^ for all ( £ 1. There exist two positive 
constants R\ and R2 such that 

Ri<Rfat)<R 2 (3.4) 

for all £ € E and for all R. 
Proof. By l|3.2|l we can write 

Ffax) = (2p+1) f ds(x + s0 2p , (3.5) 
Jo 

so that F(£,x) > for all fax) G E 2 . Moreover F(0,0) = and F(£,x) > for all £ G E if x ^ 0, 
and liiri|£|_ >00 F(£, x) = 00 for all x € E. Hence for a 7^ and 7 > 7 , one has both P(£,t) > and 
Qfa > 0, hence also R(€,t) > for all (£,t) G E 2 . Moreover lim^i^oo iZ(£,i) = 1 for all f G E, so that 
the assertion follows. ■ 



Lemma 2 Consider the equation with xo(t) the quasi-periodic solution of given in Theorem 

1. There exists 72 > such that for all 7 > 72 there is a convex set S containing the origin such that 
any trajectory starting inside S tends to the origin as time goes to infinity. One can take S such that 
dS crosses the positive and negative y- and £-axes at distances 0(7 2 ) and 0(7 2 /( p+1 )) from the origin 
respectively. 

Proof. Rescale time through the Liouville transformation 







dt' y/Rfat>),t>), (3.6) 

which is well-defined by Lemma^ Then, if we introduce the coordinate transformation -0 : (£, y) — > (v, y) 
by setting £(t) = v(r(t)) and y(t) = ■J Rfat), t) iy(r(t)), equation (|3.3() is transformed into 



v = w, 



w 



w 
R 



7 + 



R' 
2yfR 



vQ(v), 



(3.7) 



where primes denote differentiation with respect to r, Q(v(r)) — Qfat(r)) and R — i?(v(r), t(r)) 
Rfat(r)),t(r)). 

The autonomous system 

\v' = w, 
1 w' = —vQ(v), 



(3. 



can be explicitly solved: all trajectories move on the level curves of the function 

1 f v 

H{v,w) = -w 2 + dv'v'Q(v'). 

2 .In 



(3.9) 



In one has R'/VR = R/R, with 

R _ P Q 
R~ P^Q 1 

and it is easy to see (Appendix EJ that there are two 7-independent positive constants B\ and B2 such 
that 



R 



2R 



< - (B 1 +B 2 \w\ 

7 



(3.10) 
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If 7 satisfies j 2 > 2B\ we can define w as 



W = ^^>2^' (3 ' U) 

so that 7 := (7 + R'/2VR)/VR > for |io| < w. 

Consider the compact set V whose boundary dV is given by the level curve H(v,w) = w 2 /2 of the 
system (|3.8[1 . Such a curve crosses the w-axis at w = ±w — 0(^ 2 ) and the u-axis at v = 0{ r f l '^ p+1 '). 
If we take an initial datum (v(0),w(0)) G V then the dissipation coefficient 7 in (|3.7[1 . even if it changes 
with time, always remains strictly positive. Moreover H' = —^w 2 < and H' — only for w = 0, and 
for w — the vector field in 13.7J1 vanishes only at v = 0, because > for all v (cf. the proof of 

Lemma nj. Then we can apply Barbashin-Krasovsky theorem Ql^^i an d conclude that the origin is 
asymptotically stable and that V belongs to its basin of attraction. 

Let V(t) be the time-dependent preimage of V under the transformation ip. By Lemma^i! 7 is large 
enough there is a compact set S C V(t) for all fell, such that the boundary dS crosses the positive 
and negative y- and £-axes at a distances of order 7 2 and 7 2 /(p +1 ) from the origin, respectively. All 
trajectories starting from points inside S are attracted by the origin. ■ 



Lemma 3 Consider the curve g(£,t) = —7 Xo(t)) in the plane (£,?/)• There exists 73 > such 

that for 7 > 73, outside the set S defined in Lemma\^ one has 

-U 2p+1 >g(U)>--e p+1 (3.12) 

2 7 7 

for all tGt. 

Proof. Consider £ > (the case £ < can be discussed in the same way). By H3.2|) one has 



2p 

>=1 V ' / 

and, if 7 is sufficiently large so that |a;o(i)| < 2|a|, then for £ > 2|a| the sum can be bounded by 
2 2p+1 {2\a\)( 2p = C p £ 2p . Hence one has 

\e p+1 < SFfoMQ) < 2 C 2p+1 , (3-13) 

as soon as £ > 2C V (note that if £ > 2C P then one has automatically £ > 2 1 | ) . Next, we want to show 
that the latter inequality is satisfied outside S. 

Consider the intersection of the graph of g(£,,t) with dS. Let II be the quadrant {(£, y) G R 2 : £ > 
0,2/ < 0}; cf. figure H The curve dSdll is below the line 

yi(0 = yo(i-j^), (3.14) 

where y Q :— — &7 2 , with 6 > 0, is the y-coordinate of the point at which dS crosses the y-axis, and 
£0 := a^ 2 ^ p+1 \ with a > 0, is the ^-coordinate of the point at which dS crosses the £-axis. On the other 
hand the graph of <?(£, t) in II is above the curve 

2/2(0 = --(£ + 2H) 2p+1 , (3.15) 

7 

because in ^ one has £F(Z,x) < |x + £| 2p+1 + \x\ 2p+1 . 
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Figure 1: Construction used in the proof of Lemma El 



As a consequence in II the two curves dS and g(£,t) cannot cross each other for £ e (0, with 
/3 < 2/(2p +1). The latter assertion can be proved by reductio ad absurdum. First note that 07^ < £0 
for 7 large enough. Suppose that there exists £ < 07^/2 such that yi(£) = J/a(£)' Then one has 

4 < & 7 2 fl- f ) = - (e + 2|a|) 2p+1 < ^max{l2\a\} 2 ^\ 
z V 4o/ 7 7 

that is 67 s < 8 max{(, 2|a|} 2p+1 , which is not possible if j3 < 2/(2p + 1) and 7 is large enough. 

Therefore in II the graph of <?(£, t) can be outside S only for £ > aj 13 /2, which is greater than 2C P 
for 7 large enough. Hence is satisfied outside S, so that (|3.12|) follows. ■ 



Lemma 4 Consider the equation ]S.1\) . with xo(t) the quasi-periodic solution of 11. 1\) given in Theorem 
1. There exists 74 > such that for all 7 > 74, if z £ S,then either ip(t, z) enters S or crosses the y-axis 
outside S in a finite positive time. 



Proof. Consider the four quadrants 

I 

II 
III 

IV 



{(£,y)eiR 2 

{(£,y)elR 2 
{(£,?/) eK 2 



£>0,y>0}, 
£>0,y<0}, 
£<0,y<0}, 

e<o,y>o}. 



(3.16) 



In I one has £ > 0, y < 0, in II one has £ < 0, in III one has £ < 0, y > 0, and in IV one has £ > 0. It is 
easy to see that each trajectory starting in I enters II and each trajectory starting from III enters IV in 
a finite time (see Appendix |B|) . 

Consider now an initial datum z in II but not in to S. Let C\ be a continuous curve £ — > y(£) in II 
such that y < for z in II above C\\ cf. figure El Existence of such a curve follows from Lemma which 
also implies that C\ is decreasing outside S (see Appendix Q ■ The curve C\ divides II into two sets Ha 
and lib, with Ha above lib. Denote by T\ and T2 the parts of Ha and lib, respectively, outside S. Hence 
for z 6 T\ the trajectory <p(t,z) either enters S or enters T2. In the latter case it cannot come back to 
Ti, hence y{t) < y, if (£,y) = C\ n dS. This means that if the solution does not enter S then it has to 
cross the vertical axis and enter III. 

Analogously one discusses the case of initial data z in IV, outside S: their evolution leads either to S 
or to I. Hence the lemma is proved. ■ 
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Figure 2: Construction used in the proof of LemmaQ] 



Lemma 5 Consider the equation ]S.1\) . with xo(t) the quasi-periodic solution of 11. 1\) given in Theorem 
1. There exists 75 > such that for all 7 > 75 and for all z ^ S on the vertical axis, either tp(t, z) enters 
S or the trajectory ip(t,z) re-crosses the vertical axis at a point z\ which is such that \z\ — \z%\ > 5 for 
some positive z-independent constant S. 

Proof. Fix an initial datum z ^ S on the vertical axis. This means that at t = one has z — z(0) = 
(0,y(0)) outside S. Assume for concreteness y(0) > and set y(0) = l/e p+1 , with e > 0. As z(0) £ S 
there exists a constant C 2 such that 7 2 e p+1 < C 2 . Consider the change of coordinates 

X = e£, Y = e p+1 y, T = L (3.17) 



Then the system 1)3.311 becomes 



X '~ Y > (3.18) 
Y = -je p Y - XF(X, sx {e p T), 



where primes denote differentiation with respect to T. Note that Y(0) = 1 and XF(X,exo(e p T) 
X 2 P +i + 0{eX 2p ). Call S the image of S under the transformation (j3~T7|) : cf. figure 
We can rewrite the system (|3.18H as 



X' = Y, 

Y> = (X, Y) = *x(A, Y) + <F 2 (A, Y) + V 3 (X, Y), 



(3.19) 



with 

**i(X,Y) = -XF(X,sx (0)), 

* 2 (X, Y) = —X [ F{X, ex {e p T)) - F{X, ex o (0)) ] , ( 3. 20 ) 

y 3 {X,Y) = - 1£ p Y. 

If we replace ^(X, Y) with ^i(X, Y) in (|3.19() the trajectory moves on the level curve T = {(X, Y) e 
K, 2 : H(X, Y) = 1/2} for the function 

H(X,Y) = \y 2 + [ dX XF(X, ex o (0)) (3.21) 
2 Jo 
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Figure 3: Construction used in the proof of Lemma 



and crosses the vertical axis at the point (0, —1), hence at the same distance from the origin as at t = 0. 
By Lemma 13 there exists in II a curve C 2 , decreasing outside 5*, such that 

-je p Y > *(A",Y) > -~/e p Y/2, 

for Y in II below C 2 (see Appendix EJ). Such a curve can be chosen in such a way that it crosses the level 
curve r in a point P = (Xp, Yp), with Xp — 2Di(^e p ) 1 ^ 2p+1 \ for some constant D\ (see Appendix IC)) . 
Note that the time T\ necessary to reach such a point is of order 1. 

If we take into account the component ^ 2 (X, Y) of the vector field in (|3.19|) . we can move from P at 
most by a quantity of order e p+l . Indeed, as long as the motion remains close to that generated by the 
vector field ^x(X, Y), one has 

\F(X,ex (.e?Ti)) - F(X,ex (0))\ < D' 2 e \x {t?Ti) - x {Q)\ < D'{e p+ \ (3.22) 

for suitable positive constants D' 2 and D 2 ', so that the points reached at time T\ by moving according 
to the vector fields ^\ and + \&2 cannot be more distant than Z?2£ p+1 for some constant D 2 - This 
follows from the fact that the system is quasi-integrable, so that in a time of order 1 the action variable 
can change at most by a quantity of order of the perturbation as bounded in (|3.22() ; see Appendix iDl 

Finally the component \&3(X, Y) points inward along the full length of the curve T. Define T 2 as the 
time at which the the trajectory of the full system (|3.19(l crosses the curve C 2 in a point Q near P. Of 
course T 2 is near Ti, and so is of order 1, and X(T 2 ) > Xp/2 by construction, while Y(T 2 ) > Yp — D 2 e p+1 , 
with Yp > -1. 

From time T 2 onwards, we have 

I X ' = r ' (3 23) 

\Y' > ~~/e p Y/2, v ' ' 
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as long as the motion remains below C 2 . The latter property is easily checked to hold (see Appendix 
El- Then the trajectory crosses the vertical axis and meanwhile, at least, moves upward in the vertical 
direction by a quantity je p X P /2 = Dije p (-ys p ) 1 ^ 2p+1 '> . 

Therefore when the trajectory again crosses the vertical axis, this happens at a time T3 such that 
Y(T 3 ) >Y P - D 2 e p+1 + L>i7£ p (7£ p ) 1 / (2p+1 ) > -1 + Ay, with 

AY = Z?i 7 e p ( 7 £ p ) 1/(2p+1) - D 2 e p+1 > D 2 e p+1 , 

where the latter inequality holds provided ( 7 £P) 1+1 /( 2 p+ 1 ) > 2D 2 e p+1 , that is provided 

7 2(p+i) > D e p+1 , D = (2D 2 ) 2p+1 . (3.24) 

Since e p+1 j 2 < C 2 , inequality (|3~2l) is satisfied if 7 2 (p +1 ) > D C 27 ~ 2 , which requires 7 2 (p+ 2 ) > D C 2 , 
that is 

7 > (D C 2 ) 1/2(p+2} , D a = (2D 2 ) 2p+1 . (3.25) 

Under this condition one has |r(0)| - |F(T 3 )| = 1 - \Y(T 3 )\ > AY > D 2 e p+1 , so that, in terms of the 
original coordinate y, one has \y(0)\ — 1 2/(^)1 > D 2 . 

Then, if the trajectory crosses the vertical axis once more in the positive direction (and this necessarily 
happens if it does not enter 5, by Lemma 0}, this occurs at a time £4 such that 

Iv(0)|-|tf(t4)|>2£> a , (3.26) 

where we recall that the constant D 2 is independent of the initial datum y(0). Simply one can repeat the 
argument above by taking (0,y(t 3 )) as initial datum and calling £4 the time of crossing of the positive 
£-axis. This means that the trajectory either enters S or, after a complete cycle, moves closer to the 
origin by a finite positive distance 8 = 2Do- ■ 

Lemma 6 Consider the equation withxo(t) the quasi-periodic solution of J^l.l}) given as in Theorem 

1. There exists 7 6 > such that for all 7 > 76 for all z ^ S there is a finite time t{z) such that 

Proof. Consider z ^ S. By Lemma 0] either ip(t, z) enters S or there exists a time t\ such that (f(ti, z) 
is on the vertical axis outside S. Hence, without loss of generality, we can consider only initial data 
z = (£,?/) outside S such that £ = 0. Assume y > (if y < the discussion proceeds in the same way): 
we can apply Lemma 0] and we find that, as far as the trajectory does not enter <S, at each turn it gets 
closer to S by a finite quantity. Hence sooner or later it enters S. ■ 

Theorem follows from the lemmata above: it is enough to take 71 = max{7 , 72, 73, 74, 75, 76}, so 
that all lemmata apply. 



4 Uniqueness of the quasi-periodic solutions 

Let 7 be 7 > max{ 7 o,7i}. Then there exists a quasi-periodic solution xo(t) for the system <| 1 . 1 1) . and 
such a solution is a global attractor. 

In [7] we explicitly constructed a solution xo(t) with the properties stated in Theorem ^ Such a 
solution turns out to be Borel summable for d = 1. In general the solution is obtained from the formal 
series through a suitable summation procedure. Since Theorem implies that if there exists a quasi- 
periodic solution Xo (i) this has to be unique, we can conclude that for real 7 large enough there exists a 
unique quasi-periodic solution Xo{t) with the same frequency vector as the forcing /. By setting e = 1/7, 
in the complex e-plane, there is a solution x\{t) which is analytic in a domain containing a disk T> with 
centre on the real axis and tangent at the origin to the imaginary axis. For real 7 such a solution coincides 
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with xo(t) (as the latter is the only one), hence, by uniqueness of the analytic continuation, the function 
xi(t) is the only solution of (f 1 . 1|) in all the domain T>. In particular it is the only one which admits the 
formal expansion given by perturbation theory. 

We can summarise the discussion above through the following statement. 

Theorem 3 Consider the equation with f a non-zero average quasi-periodic function analytic in 

its argument and with u> satisfying the Diophantine condition 1^1. 3\) . There exists 70 > such that for all 
real 7 > 70 there is a unique quasi-periodic solution xo(t) with the same frequency vector as the forcing 
term. Such a solution describes a limit cycle in the plane (x, x) which is a global attractor. 

Note that the hypotheses made in Theorem ^ are more restrictive than those considered in [7] . In 
particular we excluded both polynomial nonlinearities and monomial nonlinearities with even degree. 

The first request aims to have a strictly positive function t), which was used in order to construct 
the positively invariant set S. We leave as an open problem to study what happens if the nonlinearity 
x 2p+1 in HI. 11 is replaced with 



For d = 1 and 7 = it is known that all motions are bounded, also replacing the constants a.,- with 
periodic functions 5,8,4. The same holds for d > 1 ^UJ- O ne could expect that the presence of friction 
tends to contract phase space toward some periodic solution (which certainly exists for 7 large enough, 
as proved in 0), but our results do not allow us to treat, in general, such a case. 

If the nonlinearity x 2p+1 is replaced with an even monomial x 2p , with p€l, then, under the further 
condition that /o > 0, there is a quasi-periodic solution xo(t): again this follows from 7 . In such a case 
xo(t) is not a global attractor, as there are unbounded solutions; cf. for example [3j for p = 1. Still one 
can prove that the solution found in |7] is unique, in the sense that it is the only attractor in a neighbour 
of the solution itself, and it is a local attractor. The same result holds, more generally, for any analytic 
g{x) in i|1.4|) such that g(co) = fo and g'(co) ^ for some Co € 1. A more formal statement is as follows. 

Theorem 4 Consider the equation \l-4\ ), with f , given by and g both analytic in their arguments, 

and with ui satisfying the Diophantine condition Assume that there exists cq € R such that 

g(co) = fo and <?'(co) 7^ 0. There exists 70 > such that for 7 > 70 there is a unique quasi-periodic 
solution Xo(t) which has the same frequency vector as f, reduces to cq in the limit 7 — > 00, and extends 
to a function analytic in a disk with center on the positive real axis and boundary tangent to the vertical 
axis at the origin. Furthermore, there exists 71 > 70 such that for 7 > 71 there is a neighbourhood B of 
the point (co,0), containing the orbit described by xo(t), with the property that all trajectories starting in 
B are attracted to the cycle described by xo(t) in the plane. 

Proof. The existence of a quasi-periodic solution xq (t) with the same frequency vector u) as the forcing 
was proved in 7j. As a byproduct of the proof, one can write xo(t) = cq + xi(ojt), with xi(rp) analytic 
in if) and of order e, if e = 1/7 (that is for e small enough one has |xi("0)| < C|e| for all tp and for a 
suitable C). Therefore we can write x(t) = xo(t) + £(i), where £(t) satisfies the differential equation 



2p+l 




a 2p +i > 0. 



(4.1) 



£+7£ + ^(£,zo(i)) = 0, 



(4.2) 



with 




(4.3) 



Then we can write Q4.2t as a system of first order differential equations, 




(4.4) 
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and define R(£,i) = F(£, xo(i)) / F(£, Co). It is easy to see that one has 



lim = l + 0(e), 

so that for e small enough one has R\ < R(£,t) < i?2, for two suitable positive constants R\ and i?2- 
Then we can rescale time and variables by setting 




r= / dty/R({(t'),1<), £(t)=v(T(t)), y(t) = ^R(at),t)w(T(t)), (4.5) 

Jo 

which transforms the system <|4.4I) into 

-vF{v,co), 

where primes denote differentiation with respect to r. 

If we neglect the friction term in l|4.6|) we obtain the autonomous system 

i V ' = W > (4 7) 

|y = -vF(v, c )) = -d x g(co) v + 0{v 2 ), 

which admits the constant of motion 

H{v,w) = ^w 2 + ^d x g(c )v 2 + O(v% d x g(c o )^0. (4.8) 

Hence the origin is a stable equilibrium point for 14.7(1 . and the level curves for H are close to ellipses 
in a neighbourhood V of the origin. It is easy to check that in V the coefficient of the friction term is 
strictly positive, for 7 large enough, because R' /2yR is (in V) less than a constant. Hence we can apply 
Barbashin-Krasovsky's theorem and conclude that the origin is asymptotically stable and V belongs to 
its basin of attraction. If we go back to the original variables (£, y) we find that V is transformed back to 
a time-dependent set V(t). But the dependence on t of V{t) is very weak (as R is close to 1 for 7 large 
enough), so that there exists a convex set 5 C for all tel. Hence any trajectory starting from S 
is attracted toward the origin. In terms of the variables (x,y), using once more that the solution Xo(t) 
is close to Co within 0(e), we can say that, for e small enough (that is for 7 large enough) there exists 
a neighbourhood B of the point (co , 0) such that it contains the cycle described by the quasi-periodic 
solution xo(t) in the plane (x,x), and any trajectory starting from B is attracted by such a cycle. 

In particular the solution xg(t) is the only quasi-periodic solution which tends to Co as e — > 0, and for 
e > small enough, say e < £0, it is the only one which admits the formal power expansion (|1.6|l . Such 
a solution was proved in jj] to be analytic in a domain T> containing the interval (0,£q), hence by the 
uniqueness of the analytic continuation, we can conclude that xo(t) is unique in all T>. ■ 

By looking at the proof of Theorem we see that it proceeds along the same lines of the proof of 
Lemma in Sectional By using the definitions of [Hj, Theorem 01 says that the orbit described by xq{£) 
is an attracting set (and an attractor) with fundamental neighbourhood B. 



5 The varactor equation, and generalisations 

The construction of the set B in the proof of Theorem^Jcan be improved of course in concrete examples. 
Here we want to consider explicitly the case of even monomials x 2p , with particular emphasis for p = 1. 
Hence we discuss explicitly the case p = 1, but the discussion can be easily extended to all p > 1. 
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So, let us fix p = 1. If we define the functions F(£,x) and R(£,i) as in the proof of Theorem^] we 



find 



R(£ t) = £ . + 2x o( t ) = x , r, x o(t) - Co 



(5.1) 
(5.2) 



£ + 2c £ + 2c 

where Co = \ffo, with / > 0. For 7 large enough, that is for e small enough, one has 

\xo(t) - c | < Ci£, 

for a suitable constant C\ (see the first lines of the proof of Theorem HJ. 

The transformation of coordinates described in the proof of Theorem 0] leads to (|4.6H . with 

v F(v, c ) = v 2 + 2c v, (5.3) 



so that, if we set 7 = (7 + R'/2-\/R)/yR, we can interpret the system as a Hamiltonian system with 
hamiltonian 



H(v,w) = ^w 2 + U(v), 



U(v) = ^v 3 + c Q v 2 , 



(5.4) 



in the presence of a friction term with non-constant dissipation coefficient 7. If we neglect the friction 
term, that is if we put 7 = 0, then the system becomes the Hamiltonian system 



w 



2c w, 



(5.5) 



which admits a stable equilibrium point P\ = (0,0), an unstable equilibrium point P2 = (— 2co,0) and 
the separatrix T of equation 



= ±V2(ff(-2co,0)-tfH), JT(-2q,,0) = ^, 



(5.6) 



which contains a homiclinic orbit to the right of P2. Consider the open set T>\ containing the point P\ 
and with boundary given by (the closure of) the homoclinic orbit in T. Any level curve with energy 
V(0) < E < V{— 2co) contains a bounded connected component, internal to T>±, which is a closed 
orbit Ce for the system 1)5. 5p . Call C the closed curve Ce such that it intersects the negative w-axis at 
v = — 2cq + C2s" , for suitable constants Ci and /3 to be fixed, and T> the bounded open set with boundary 
C. 

From l|5.1fl . one has 

(5.7) 



so that R(£,,t) — > 1 as e — > 0, provided /3 < 1; for instance one can take /3 = 1/2. Then, if 7 is large 
enough, for any point (v, w) &T> one has 



\R'&t)\< 



2£'(co - x ) 


+ 


2x' 


« + 2c ) 2 




C + 2c 



for suitable constants C3 and C4, so that, by using also that £ = u and 

|C + 2c | ^Cae' 3 

for all (f,w) € T), one obtains 

i?' 



17-7 



< 



2R 



< C^e 1 ^ < C 5 s 2 ^h < 7 



(5.8) 



(5.9) 



(5.10) 



for a suitable constant C5 , hence 7' > for 7 large enough. Therefore we can conclude that T> is positively 
invariant, and we can apply once more Barbashin-Krasovsky's theorem to conclude that the equilibrium 
point Pi is asymptotically stable and the set D is contained inside its basin of attraction. 
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Figure 4: Estimate of the basin of attraction (grey set) of the limit cycle (solid line inside the grey set) for large 
7. The curve To is the separatrix of the system described by the differential equations x + x' 2 = cq. 

If we go back to the original coordinates, we see that, for 7 large enough, there is a quasi-periodic 
solution xo(t) which moves very close to the point (co, 0) (within a distance of order e — I/7); cf. figure^] 
Its basin of attraction contains a large set whose boundary is at a distance of order e 13 from the separatrix 
To of the system, that is from the curve of equation y = ±-\/2(2cq — x 3 + 3cqx)/3. 

6 Improved estimates for the basins of attraction 

We start by constructing an invariant set A, valid for all p £ M, which has the property that its area is 
0(7). The construction is surely not optimal but is included because it is, amongst many sets constructed, 
one which grows with 7 in the vertical direction. Constructing invariant sets containing the attracting 
orbit with 7 considered to be large is not difficult, but finding ones that grow with 7 was found to be less 
trivial. The construction follows very closely the ideas exploited in 0. 
Rewrite (|1.4|) . with g(x) = x 2p as 

y = /M) - iv - % p , 

so that the vector field generated by the differential equation is defined by <p(t) = (y, f(u3t) — "fy — x 2p ). 
Let f 2p < f(ut) < F 2p . Note that in this and the following section, / is just a constant and not the 
function t — > f(u>t). We adopt this notation here to conform to that used in the analysis in 0. For the 
same reason we denote in boldface vectors in K 2 , so that a ■ b denotes the scalar product in H 2 . 

The two vector fields, <p F = (y, F 2p — 72/ — x 2p ) and <pj = (y, f 2p — jy — x 2p ), have no explicit time- 
dependence and also have the property that for all t, <p(t) = ^<p F + (l—^i)(pj where /i is a (time-dependent) 
scalar. 

Following , we let the boundary of invariant set A be a hexagon GHIJKL whose edges are straight 
lines, except for HI and KL. GH and JK are horizontal and LG and IJ are vertical — see figure The 
co-ordinates of points H and K are (0,yii) and (0,yi<) respectively. The dotted curves in figure El are 
Pf ■ y = (F 2p — x 2p )/'y (upper) and Pf : y = (f 2p — x 2p )/j (lower). As shown in [3], only in the region 
between these curves is the sign of y ambiguous; above Pp, y < and below Pf, y > 0. Provided that 
y H > F 2p /-f, which will turn out to be automatically satisfied, <j>{t) will be into A along GH, and point J 
being below Pf guarantees that 4>(t) will always be into KJ, both results holding for all t. Furthermore, 
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the sign of x — y guarantees that <fi(t) will be into LG and IJ, again for all t. It therefore remains to 
prove that <p(t) is into the curved sides HI and KL for all t. 



G H 




Figure 5: The absorbing set GHIJKL, A. Also shown as dotted lines are the curves Pf ■ y = (F 2p — x 2p )/ , y and 
P } : y = (f 2p - x^)H- 

Let us define HI by y — Ai (F 2v — x 2p ) with Ai <E E, positive; then the inward-pointing normal, rt\ = 
(— 2pAi£ 2p_1 , —1). In order to prove that <p(t) is into HI for all time, we need only show that n\ ■ 4> F > 
for x G [0, F]. We have ni -4> F = [F 2 p - x 2p ] [Ai( 7 - 2p\ 1 x 2 P~ 1 ) - l] and, since F 2 p - x 2 p > 0, and the 
second bracket in the scalar product reaches its minimum value over [0, F] at x — F, we require 

2pF 2p - 1 X 2 - 7A1 + 1 < 0, (6.2) 

in order for the scalar product to be non-negative. Solving this quadratic and choosing the larger solution 
gives 

7 



Ai = 



4 P F 2 p- 




8pF 2 P~ 



(6.3) 



provided Ai € M. This is true for j 2 > 8pF 2p 1 . We now show that this condition on 7 also forces 
yn = X\F 2p > F 2p /j, or Ai > I/7, to be satisfied. Letting q = ApF^ 1 , the condition for Ai to be real 
is 7 2 /<? > 2; using this in the definition of Ai gives (rf/q) (l + y/l - 2q/j^j > 2 ^1 + y/l - 2q/j 2 ^ > 1, 
which is clearly true. Hence, boundary HI has been constructed in such a way that <p(t) is into it for all 
time, provided that H is above Pf'- this is true if Ai > 1/7. Note that Ai ~ 7 / (2pF 2p ~ 1 ) as 7 — > 00, and 
so y H = XiF 2 p = 0(7). 

We now define KL by y = A2 (/ 2p — (x + 2/) 2p ) where A2 € R > 0. This has inward normal 
n,2 = (2p\2(x + 2/) 2p_1 , l) and the inequality we need to consider is now 

4> t ■ n 2 = A 2 [(x + 2ffP - f 2p ] [7 - 2p\ 2 (x + 2/fP- 1 ] + - x 2 p > (6.4) 

for x € [— /, xk] where xk, satisfying F > xk > — /, is to be defined. The first term consists of a product 
of two terms, the first of which is positive for x G (— /, 00) and the second of which is a monotonically 
decreasing function of x for x G R \ {—2/}; it is zero at x = xq = (7 /2p\2) 1 ^ 2p ~ 1 * > — 2f. Choosing 
xk = x — ensures that the product term in (|6.4|l is non-negative and gives 

A2= 2 P (2/V r (6 ' 5) 
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The last term in inequality (|6.3|l is non-negative for x € [— /, /] and so the scalar product is non-negative 
for A2 as given above, provided that xk = 0. We also require point J to be below pf. The y-co-ordinate 
of J = yK = A2 (/ 2p — (2/) 2p ) and so this condition becomes 



^ ^ p (F 2p - f 2p ) 



/(1-4- 



p) 



(6.6) 



We can now state our result. 



Theorem 5 Define A as the hexagonal set whose vertices are G, H, I, J, K and L, where 

G=(-f,\iF 2p ), H=(Q,\ 1 F a *), I=(F,0), 
J=(F,-X 2 f 2 P(4P -1)), if= (0,-A 2 / 2p (4 p -l)), L=(-f,0), 



(6.7) 



and where Ai and A2 are given by equations \b\ and IB^I), respectively. Let the edges LG, GH, IJ and 
JK of A be straight lines and let HI be given by y = Ai (F*p~— x 2p ) and KL, by y = A 2 (f 2p - (x + 2/) 2p ) . 
Then, provided that 

„_! p{F^~pP)\ 



7 2 >max^ 2p - , /(1 _ 4 _ p) j, 

set A is an invariant set containing the limit cycle described by the solution xo(t) in the plane. 
Figure compares the set GHIJKL with the actual basin of attraction for p — 1, 2 and 3. 

80 



(6.8) 



PS 




20 40 60 80 



Figure 6: The set GHIJKL compared to the set B c (cf. 0) with 1 < f(t) = (5 + 3sint)/2 < 4, 7 = 9 and p = 1 
(left), p = 2 (middle) and p = 3 (right). The set £? is defined such that all trajectories starting outside B are 
(numerically) found to remain outside, hence it provides an estimate for the basin of attraction of the periodic 
solution, xo(t). 

The existence of the set A, with the properties stated in Theorem El allows us to improve the estimate 
T> of the basin of attraction found in Section for the varactor equation studied in , that is for the 
system (|7.1H . where f(u)t) = a + /3sint, with a > and \/3\ < a. Define Vq = VUA. First note that both 
D and A are strictly contained in 2?0i a s follows easily from the remark that T>q is inside the separatrix 
Fo (cf. the end of Sectional), and at a distance 0(1/^/7) from it. The inclusion D C T>q is obvious 
as A grows linearly at 7 in the vertical direction. The inclusion A % V follows from the fact that the 
separatrix intersects the ir-axis in x = — cq and x = 2co, with cq = -y/a; hence xl = —\fa — p > -co and 
xi = y/a + {3 < \plct < 2co- Then the set T>i := ^4 \ (.4 n T>) is non-empty. Furthermore a trajectory 
starting in T>\ can go out of T>\ only by entering T>q (by the invariance of A). This means that all 
trajectories starting in T>q are attracted by the limit cycle described by Xo(t). 
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7 Blow up in finite time 



We prove finite-time blow up for the system l|6.1|l by first finding an invariant set J (Lemma 0) and then 
constructing an invariant subset, S(— Xo), of J ', the latter construction being carried out in Lemma El 
Within S(—Xo), a differential inequality must hold and we show in Theorem that all solutions of this 
differential inequality must blow up in finite time. 

Lemma 7 Define h(x) = 2px 2p - 1 (F 2 p - x 2 p) - 7 2 (F 2 p - f 2 P) and the set 

V = {(x, y)\x < {F 2 p - x 2 P)h < y < 0}, (7.1) 

where — £ is the root of h(x) with — £ < — F. Then £ as defined exists uniquely and T> is an invariant set. 

Proof. The set D is clearly absorbing along its horizontal, y = 0, x < — £, and vertical, x = — £, 
(^2p _ £ 2 p)/j < y < boundaries. We therefore only need to prove that it is also absorbing along the 
curved boundary, y = {F 2p — x 2p )/ 7 for x < — £. The appropriate normal here is n = (2 P x 2p ' 1 /j, 1) and 
so we require n ■ cf>f > 0, which is equivalent to proving that h(x) > 0, for x < — £. 

To this end, first note some elementary properties of h(x): (i) h(—F) = 7 2 (F 2v — / 2p ) < 0; (ii) h(x) 
has exactly three stationary points for x € R: /i'(ir) = for a; = and a; = x± = ±i^[(2p— l)/(4p — l)] 1//2p ; 
(hi) h"(x-) > 0; and (iv) h(x) — * oo as x — > — oo. From these, it becomes clear that — F < x_ < and 
ft.(x_) is a minimum; and, in the light of (i) and (iv), h(x) has exactly one real root, — £, G (— oo, — F). 
Hence, /i(x) > for x < — £ and the invariance of set V is proved. ■ 

We now define, for x < — Xo, a curve G(x, y) : y = —b (—X — x) p with — X < — £ < 0, p = 3/2, and 
b > to be found. In order that <S(— Xq) C J7, — Xo < — £ and > 0. In the proof of Theorem El it 
will be required that S(—Xo) has infinite area for all p G IN, and hence p < 2 — otherwise, the curves 
(F p — £ 2p )/7 and G could intersect at some finite x < — £ when p = 1. In order that solutions blow up 
in finite time, p > 1 will also be required: hence, we choose p = 3/2. 

The proof of the invariance of S now follows. 



Lemma 8 The set 

where 



S(-X ) = {(x, y)\x < -Z, U > y > (F 2p - x 2p )/ 7 }, (7.2) 

U= \ -b(-X ~x) 3/2 x<-X ^ 3 % 

1 otherwise 



with — Xo < — £ and 6 sufficiently small, is an invariant subset of J , and contains points (x, y) with 
y -> -oo. 

Proof. The fact that the vector field is into all the boundaries of S except y = —6 (— Xo — x) 3 ^ 2 has 
been proved in Lemma therefore, we only need to consider the boundary y = — b{— Xo — x) 3 ^ 2 . The 
appropriate normal is that which points into S, this being n = (3b/2)y/— Xq — x, —1). The correct choice 
for the vector field here is <p F = (y, F 2p — jy — x 2p ), which gives 

n ■ <t> F = -b(-X - xf/ 2 [(36/2) V-X -x + 7] + x 2p - F 2p > 0, x < -X . (7.4) 

Since — Xq — x > 0, we can substitute v 2 = — Xq — x, giving 

L(v) = (Xq + v 2 ) 2p - F 2p -bv 3 (3bv/2 + -f) > 0, v > 0. (7.5) 

Let l(x) — x p ; then the mean value theorem states that, for b > a, 1(b) — 1(a) = (b — a)pc p ^ 1 , where 
c £ (a, b). Applying this to the first two terms in equation Q7.5JI . we have 

L(v) > [(X Q +v 2 ) 2 - F 2 ]pF 2p - 2 -bv 3 (3bv/2 + ~f) 

= v\ V F 2 p- 2 - 3b 2 12) - 7 bv 3 + 2X oP F 2 P- 2 v 2 + (X 2 - F 2 )pF 2 P~ 2 

> v 2 [v 2 (pF 2p - 2 - 3b 2 /2) - jbv + 2X Q pF 2p - 2 ] = v 2 M(v), (7.6) 
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where we have used Xq > F to obtain the last inequality. Hence, we need to show that for b > 
sufficiently small, M(v) > for v > 0. First, write M(v) = aiv 2 — a\bv + ao, ai,aa > 0, and let b 
be sufficiently small that 02 > 0. Let M(v) = M(v)/a.2 = v 2 — 2c\bv + cq, with c\,cq > 0. Then 
M(v) = (v — bci) 2 + co — bcf, and it is plain that b can be chosen to be small enough that M(v), and 
so L(v), are non-negative for v > 0. 

We will also require <S(— Xq) to extend to infinite negative y values; it has this property provided that 
the curves y = —b (—Xq — x) 3 ^ 2 and y — (F 2p — x 2p )/^ nowhere intersect for y < 0. This, too, is clearly 
true if b is small enough. ■ 

It is possible that some of the conditions applied in the above proof could be relaxed, but a 'better' 
invariant set S is not required in the proof of the following theorem. 

Theorem 6 For all Xq = (xq^q) G J , 3 — Xq < — £ which is such that Xq G S(—Xo). All solutions 
starting from such an Xq blow up in finite time. 

Proof. Let Xq — (x(0),y(Q)) = (xo,yo) G 3 ■ Then there always exists an — Xq < — £ such that Xq is in a 
subset S(— Xq) of J: choose any — X > x . Also, since x < in J, x < xq and so — X — x > — X — x > 
0. Additionally, since x a G S(—Xq), the differential inequality — b(— Xq — x) 3 ^ 2 > y = x > (F 2p — x 2p )/~f 
applies for all time t > 0, by the invariance of S(— Xq). 

In fact, only the upper bound is important here, and with the substitution u = —Xq — x > and 
uq = —Xq — xq > 0, this becomes u > bu 3 / 2 . Integrating gives 

dipip- 3/2 >b f ds, (7.7) 

u Jo 



-1/2 



— u 1 / 2 > bt/2, which, after re-arranging, gives 



u > 



u 1/2 -bt/2 for t G [0, too), (7.8) 



— 1/2 

where btoo = 2u . The above inequality shows that x(t) tends to — oo within a finite time t < t^. By 
the invariance of S(—Xq), x(t) cannot tend to — oo without also y(t) — * — oo, and so finite time blow up 
is proven. ■ 

Note that for p = 1 there is a set B (cf . figure [fjj such that all solutions starting from B reach J in 
a finite time: this was proved in Therefore for p = 1 Theorem [H] shows that all trajectories starting 
from B blow up in finite time. 



8 Conclusions, extensions and open problems 

We conclude this paper with a list of open problems (some already mentioned in the previous sections). 

The first one concerns possible extensions of the proof of Theorem [3] to the case of more general 
polynomials of the form 1)4.1(1 . A natural question is under which conditions there is still a global 
attractor, in these cases, when the dissipation coefficient is large enough. 

A characterisation of the set B can be given in some concrete cases, such as that of the varactor 
equation considered in Section [5J Its diameter is of order 1. We have seen in Section that we can 
improve the estimate by obtaining a set whose size increases linearly in 7 in the vertical direction, but in 
such a way that it is still expected to be strictly included inside the actual basin of attraction. It would 
be worthwhile to attempt constructions of sets contained inside the basins of attraction that are as large 
as possible. 

We also leave as an open problem for the varactor equation the proof that any bounded solution is 
attracted by xo(t). On the basis of numerical simulations, we conjecture that this is the case. 
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Another interesting problem is whether one can weaken the hypotheses on the function g, both for 
determining the existence of a quasi-periodic solution with the same frequency vector as the forcing and, 
in that case, for proving its uniqueness and attractivity. 

Finally, extensions to higher dimensional cases would be desirable. 

A Proof of flHUnD 

One has P = d^P £ + d t P and Q = d^Q £, so that 

P Q_(d^P d tL Q\ ■ | d t P djP djQ | djQ {Q | d t P 



P Q \ P Q ) P P PQ ' P 

One can write 

2p , 

d * P = E( t l^-Vo-^tJioW, (A.2) 

3 = 1 V 

'2p+l 



3 = 1 
2p 



Q-P = £( 7 ^^(a-^oW) 



3 = 1 

where Xo(t) — 0(1/7) and xo(t) — a — 0(1/7). 

Finally |x (f)| < 2\a\ for all f G E if 7 is large enough, and both ^P+^/p and £ 2p+1 - j /Q tend to 
zero as £ — ► 00 for j > 1. Hence ()3.10|) follows, with the constants B\ and -B2 depending on p but not on 
7- 

B Initial data in I and III 

Take an initial datum z = (£, y) in I. If y = then £ = and y = — £F(£, xo(t)) < 0, so that the trajectory 
enters II. If y > then £ > and y < 0. 

Moreover 9 s (£F(£,a;)) = <9 ? (x + £) 2p+1 = (2p + 1) (x + £) 2p > for all x G E, so that, by using the 
fact that £(£) > £(0) as long as (£(f), y(i)) remains in I, one has in I 

ZF(Z,x {t)) > wiZF(S,x (t)) > m infF(Z{0),x {t)) > c> 0, (B.l) 

where we used that F(£, x) is strictly greater than a positive constant for x ^ (see the proof of Lemma^). 
Therefore we obtain 

y<-iv- c, (B.2) 
which implies that y(t) reaches the £-axis in a finite time. 
Analogously one discusses the case of initial data z in III. 

C On the curves C\ and C2 

Call T the subset of II outside S. 
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Define C± as a continuous curve in II such that in T it is given by the graph of the function £ — > 
— £ 2p+1 /4j. In (|3.3|) one can write y = j(—y + <?(£, t)), with <?(£, t) defined in LemmaEl By Lemnia|3 in 
T one has g(£, t) < —^ 2p+1 /2j, so that at all points in T above C\ one has 

-v+g&t)) = \v\+ g &t) < -±e^ < -±e»\ 

hence y < 0. 

Define C 2 as a continuous curve in II such that in T it is given by the graph of the function £ — > 
-4£ 2p+1 /7. By Lemma El one has i) > -2£ 2p+1 /7, so that in all points of T below C 2 one has 
j/ < -4£ 2 p +1 / 7 < 2p(^,t), hence 

-TV > i{-y + 9(€, *)) > -7J//2, 
so that —72/ >y> —jy/2. In terms of the rescaled variables (X, F) this yields y' = $f(X, Y), with 

_ 7£ py > «jr(X, y) > - 7 e p y/2, (C.l) 

as asserted after (|3.21|l . 

The point P is given by the intersection of the curve Y\{X) = —4X 2p+1 / r ye p withe the level curve T. 
Hence 

1 1 / AX 2p+1 \ 2 nip+i 

2 = 2 {^y) + XF (X,e Xo (0)) = + 

hence X P = 0{{- i e p ) l / {2p+l ^). 

Now consider the solution of (|3~Tg|) with initial datum Z(T 2 ) = (X(T 2 ),Y(T 2 )). We want to check 
that the solution remains below C 2 until it crosses the y-axis. The solution of 

| X ' = r ' (C2) 
\Y' = V{X,Y), V ' ' 

with &(X, y) satisfying the bounds (jC.lf) . moves below the line with slope — je p passing through Z(T 2 ), 
that is below the line of equation 

Y 1 {X)=Y Q - 1 e p X (C.3) 
with y determined by the request that for X = X = X(T 2 ) one has 

Y - 7 s p X = -^-X 2p+1 , (C.4) 

7£ p 

where the graph of —4X 2p+1 /je p describes the curve C 2 in the coordinates (X,Y). By using that X 
is close to Xp one realises that Yq in l|C.4(l has to be negative. In turn this implies that the line of 
equation (|C.3J| is below the curve C 2l so that also the assertion after H3.23J) is proved. 

D Variations in finite times for quasi-integrable systems 

The system obtained from l|3.19(l by replacing ty(X,Y) with ^>\(X,Y) is an integrable Hamiltonian 
system, with Hamiltonian (|3.21(l . For e — the Hamiltonian reduces to 

H (X,Y) = l -Y 2 + ^-^X 2p + 2 , (D.l) 

which can be written in terms of the action-angle variables (J, ip) as 

H a (X,Y) = H (I) = Cp I {2n+2)/{n+2 \ (D.2) 
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where c p is a suitable p-dependent positive constant. By taking into account the other terms of the vector 
field, we obtain 

H{X, Y) = H(I) = Cp /(2n+2)/(n+2) + 0(/ (2„+l)/(„+2) ) _ ^.3) 

The equations obtained by adding to ^i(X, Y) the vector field ^(A, Y) are still Hamiltonian, and are 
described by the non-autonomous Hamiltonian 

H{I)+H x {I,v,t), (D.4) 

with Ha given as in l|D.3(l and Hi of order e p+1 as long as the action variables remain of order 1. 
The corresponding equations of motion are 

i = -d ¥ H l {i,ip), 

tp = uj (I) + dMiil,^), 

with lvo{I) = diHo(I). Then one immediately realises that in a time of order 1 the action variables 
remain close to their initial values. In turn this implies that also the angle variables are changed by order 
e' p+1 with respect their unperturbed values. In terms of the original coordinates (X, Y) this means that 
the solution remains within a distance 0{e p+1 ) with respect the unperturbed value. 
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